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The notion of self-similar energy cascades and multifractality has long since been con- 
nected with fully developed, homogeneous and isotropic turbulence. We introduce a 
number of amendments to the standard methods for analysing the multifractal proper- 
ties of the energy dissipation field of a turbulent flow. We conjecture that the scaling 
assumption for the moments of the energy dissipation rate is valid within the transition 
range to dissipation introduced by Castaing et al.{ Physica D 46, 177 (1990)). The 
multifractal spectral functions appear to be universal well within the error margins and 
exhibit some as yet undiscussed features. Furthermore, this universality is also present 
in the neither homogeneous nor isotropic flows in the wake very close to a cylinder or 
the off-centre region of a free jet. 



1 Introduction 

The concept of an iterated distributive process behind turbulence, commonly called 
the "energy cascade", is much older than fractal geometry. It might be said that 
multifractal theory has its roots, partially, in turbulence research Id. Of course, by 
being finite the statistics of the energy cascade do not have to be multifractal, and 
even if it may be conceived to be so, the expedience of a multifractal description 
is far from obvious. For several decades one incentive for turbulence research has 
been the idea of a universal behaviour in small scale intermittency 13, but it has not 
been possible to bring this universality to light in the regime of length scales which 
is accessible to a multifractal model— JChese questions were adressed in pioneering 



works by Meneveau and SreenivasanH'Q'O, and it became clear then that multifractal 
geometry is an issue in the long process of understanding turbulence. 

This paper is devoted to an extensive examination of experimental data in a 
multifractal framework. We shall not venture to relate our experimental findings to 
multifractal models of turbulence here, however, we vouch for the multifractality of 
turbulence. 

The central quantity of our investigation, as introduced by Obukhov i and 
Kolmogorov Q in 1962, is the energy dissipation averaged over a volume of size r. In 
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the case of locally homogeneous and isotropic turbulence this is approximated by 

where u{x) denotes the longitudinal velocity component at point x along a 1-D cut 
through the energy dissipation field (Figure 0) . 

We define the multifractal measure fii on a grid of size r to be fii — rer{ir). 
The q-th order moment of /i is defined as : 

i 

where the sum is taken over all boxes of some r-grid and T{q) is the multifractal 
scaling function. 

To derive these quantities from the data, we employ an improved multifractal 
box counting algorithm recently suggested eI. We sketch this algorithm in section 
two. Section three concerns the details of the experimental data and the particu- 
lar aspects of the selection of the proper scaling range. Section four presents the 
results obtained for locally isotropic flow conditions indicating a universal multi- 
fractal behaviour in the range defined by the transition from the inertial range to 
the dissipative limit. These universal features are also recovered in inhomogeneous 
and not fully developed turbulence, as shown in the concluding section. 
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Figure 1: The energy dissipation time series (equivalent to a 1-D cut along the mean flow) for 
experiment Jl. Units on the x-axis are in terms of the Kolmogorov microscale r). 



2 An Improved Multifractal Box-Counting Algorithm 

The necessity of taking great care when utilising fixed-size box counting algorithms 
has frequently been expressed. However, here we deem fixed-size algorithms to 
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teethe method of choice for two reasorLai—thcir shortcomings are well understood 
□'Ej and can in some cases be balanced ISO, and they are dnser to the theoretical 
concepts of turbulence than their fixed-mass counterparts 113 in that they do not 
mix length scales . We also refrain from the direct evaluation of the /IjjiL curve 
E3, because we find artefacts are easier to trace in the standard approach 0O. The 
information content of both methods is equivalent since they are based on essentially 
the same counting procedure. 

We want^to outline in brief some of the improvements introduced to the method 
of moments □. The prime source of errors was identified as the finite size artefacts 
due to ill-fitting grids. To reduce the influence of these "clipping errors" we replace 
in eq. (H) the local measure fii by a local average < /i >i and define 

n 

< f^>;^ ^^^'[[Q{^^{X^+SJ)) for q>0 (3) 

i=i 

</u>»= H'^^^^ + '^j) 9<0 , (4) 

where {6j} is a set of displacements smaller than r. 

Secondly, we introduce the concept of wandering intercents-which describe the 
deviations from linear scaling of the partition function (|^) B'tSEj. For a given q we 
define 

C{q,r)^ log S,{r)-T{q)\ogr , (5) 

where f{q) is an estimate of T{q), e.g. the result of a least squares line fit to 
logS'q(logr). The C{q,r) can be thought of as a "lacunarity jfunction" . Note that 
the concept of Extended Self-Similarity introduced by Benzi Ell is a special case of 
eq. (H). Our development relies on the assumption that the C{q, r) does not behave 
too erratically in q and can ideally be factorised into functions of r and q only. 
Although we shall not make detailed assumptions about the C{q,r) in general, a 
lack of coherence with respect to g at a noint go betrays the presence of a phase 
transitional behaviour of the multifractal B. For these purposes we employ plots of 
C{qi,ri) against C{q2,ri) for qi < q^ < q2 and expect that a lack of correlation 
between the fluctuations in r of the wandering intercepts appears as a scattered 
plot. 

The greatest obstacle in the multifractal analysis of turbulence data (as with 
any other physical data) is the determination of the proper scaling range. For our 
findings it is characteristic that the multifractal scaling assumption applies only to 
the transitional range to dissipation, which can loosely be defined by the lengthscale 
r] where the local Reynolds number becomes of the order of unity: ^ — I. In most 
cases, the lower bound is dictated by the detector resolution, which usually exceeds 
77 quite considerably for turbulent flows of a reasonable Reynolds number. The 
upper bound on the scaling range is much more indistinct, since typically one can 
only discern a sweeping transition to a power law behaviour of Sq (r) , which in turn 
sometimes extends over less than one order of magnitude. Although it is tempting 
to expand the scaling range both to acquire a better statistical basis and to extend 
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the validity of the analysis, we handle the problem quite restrictively. Since the 
essence of numerical multifractal analysis is extracting to a high precision minor 
differences in interrelated scaling laws, we deem our approach justified. 

3 The Statistical Properties of Experimental Data 

The turbulence velocity data were obtained in three series of experiments of flows 
in air. Series J was a free jet through a nozzle of d = 8 mm diameter and an output 
velocity of 40 m/s. The estimated R\ in the centre was 210. The jet was scanned 
axially and radially at an axial distance of 60d = 48 cm. The time series had a size 
of 1.25 X 10'' data points. Series CF was a wind tunnel experiment with a cylinder of 
diameter D = 5 cm. The turbulent wake was examined at a distance of 32d = 160 
cm for a number of velocities, yielding R\ from 200 to 520. The size of the time 
series was 1.25 x 10^. Series CV consisted of a variation of the location of the 
probe in the wake field of the cylinder at 21m/s velocity, ranging from Ad = 20cm 
to AOd = 200cm. Again, 1.25 x 10^ points were sampled. The equipment was a 
DANTEC anemometer StreamLine 90CN10 with single wire probes 55P01 (sensitive 
length 1.25mm) for the jet measurements and x-wire probes 55P61 (sensitive length 
1mm) for the cylinder measurements. The use of the latter probe facilitated a more 
accurate measurement of the longitudinal velocity component. 

An extensive analysis was performed for all time series. This included the 
calculation of characteristic mean-field quantities (see table 1) and the generation 
of probability density functions (PDFs) of velocity increments. The changing form 
of the PDFs as a function of ihe spatial resolution r was quantified by a parameter 
A^(r) according to Castaing iB. The power law behaviour of the resulting function 
A^(r) locates the_transitional range between the classical inertial range and the 
dissipative limit t3 (see figure |^). In the following we restrict our investigations 
of the scaling behaviour of < >, respectively Sq{r), to this transitional range. 
Critical evidence for the coincidence of these scaling ranges was deduced from our 
experimental data. The chosen ranges for each example are given in tables 1 and 2. 

It is one of the traits of box-counting, that it delivers a result under almost 
all circumstances. However, it is sometimes open to arbitrariness to find out which 
results comprise significant information about the multifractal. As it turned out, 
statistical resolution is a severe issue for moments as small as 3 for about 1 million 
data points, |_and 10 million points improve the situation only marginally. It has 
been shown EJ, that /(a) will assume negative values due to the 1-D nature of the 
data set as opposed to the three dimensional energy field, which infers that the 
main contribution to a moment that corresponds to a negative /(a) comes from 
very few singularities. This being the case, it comes as no surprise to find that 
the correlation plots for qi,q2 > 5 are extremely well correlated. In general, we 
attribute little credibility to any moment larger than 4. For negative q, noise and 
digitalisation artefacts put an end to the reliability of the results at sometimes quite 
moderate values oi q ~ —1. A special case are the inhomogeneous examples, where 
due to the background of the intertwined laminar flow no multifractal fine structure 
can be resolved with negative moments. 
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Figure 2: A plot o/log against log(r/r;) for experiment Jl. The transitional range corresponds 
to the linear part of the curve abutting the maximum. 



4 Multifractal Spectra in the Locally Isotropic Case 

In the following, four examples of classical fully developed turbulence are inves- 
tigated in the centre of the flows, where local isotropy is known to hold for small 
scales. The experiments were chosen to cover a range of Rx from 190 to 520. Table 1 
gives the relevant details. An example of an energy time series is given in figure |l|. 

In the context of this paper, we want to establish experiment Jl as the standard 
for a simple reason: the detector resolution and rj are quite close, so that one can 
assume that the fine structure of the original velocity signal is well resolved. In 
fact, a direct comparison with experiment J2 shows that this point is crucial for the 



calculation of extreme moments. Figures |3.1| and |3.2| show the T{q) and f{a) curve 
for the transitional range. 

Figure ^ shows a superposition of the respective /(a) spectra for all experi- 
ments of table 1. The errorbars are significantly larger than the deviations between 
the curves. Given these results, we have to assume universality of the spectra of 
singularities. 

Although it is notoriously difficult to establish negative values of /-reliably, we 
have confidence in the results we have given. A self consistent analysis El of the T(g) 
curve showed that statistics are still sound for q ~ 3, whereas the intersection with 
the a-axis occurs for q = 3.2 c.f. figure ^ 

There is a couple of interesting, but less obvious results. Firstly, from correlation 
plots of C{q,r) one is led to deduce the presence of a phase transition at go — 1, 
which does not express itself in a kink in /(a) since both tangents necessarily have 
slope 1. We interpret this as the presence of two phases: the highly intermittent or 
turbulent phase which dominates the moments for q > 1, and the "laminar" phase 
which does so by the less pronounced singularities for moments with q < 1. 

Secondly, if one sets a threshold for the analysis of the measure, one finds 
that for all (very small) thresholds the laminar branch of the multifractal spectrum 
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3.1: The Experiment Jl serves as a stan- 
dard against which the others can be com- 
pared since the detector resolution is close 
to rj. As it turns out, resolution has an im- 
pact both on positive and negative moments. 
This figure shows the T{q) curve of this ex- 
periment 



3.2: The /(a) curve cor resp onding to exper- 



iment Jl, see also figure 3.1. This curve was 
obtained by a Legendre transform of T{q). 
Errors are amplified due to the great lever- 
age of q. 



disappears. The conclusion is, that aU energy dissipation must be concentrated 
on the set which corresponds to the Holder exponent a{q = 1). Whilst this is 
necessarily so for mathematical multifractals, it is remarkable for a physical, hence 
finite one. In other words, the steepest singularities feed on most of the energy and 
dissipation attracts further energy to dissipate. 



flow 


free jet 


free jet 


wake 
of a cylinder 


wake 
of a cylinder 


reference 


Jl 


J2 


CFl 


CF2 


position of probe 


125 diam. 


60 diam. 


32 diam. 


32 diam. 


number of points 


1.25 X 10^ 


1.25 X 10^ 


1.25 X l{f 


1.25 X 10^ 


mean velocity 


2.26 m/s 


4.2 m/s 


22.2 m/s 


4.6 m/s 


sampling rate 


8 kHz 


50 kHz 


120 kHz 


20 kHz 


Kolmogorov 
microscale rj 


0.25 mm 


0.1 mm 


0.07 mm 


0.22 mm 


Taylor microscale A 

> _ {{u-uf) 
^ - {(dAu-u))-) 


27 ?7 


28 ri 


44 T] 


29 77 


Taylor-Reynolds number 


190 


210 


520 


220 


transitional range 


4 .. 113 -q 


8 .. 83 


26 .. 79 r; 


8 .. 85 r; 


detector resolution 


4 rj 


10 Jl 


14 


5 T] 


permissible moments 


-2 .. 5 


-0.2 .. 4 


-0.6 .. 4 


-4 .. 4 



Table 1 Summary of experimental conditions and evaluation parameters for the 
locally isotropic flows. 



5 Multifractal Spectra in the Inhomogeneous Case 

Inherent to most theories pertinent to small scale turbulence is the assumption of 
isotropy and homogeneity of the turbulent velocity field. A question of practical 
importance is if similar properties can be found where these conditions are violated. 
To approach this we investigated a free jet with an radial offset of 9 cm in 48 cm axial 
distance from the nozzle, i.e. in the outer layers where laminar phases are mixed 
in. As a second inhomogeneous situation the wake of a cylinder was examined 
in 4 diameters distance, i.e. under conditions of not fully developed turbulence. 
Figure || shows the energy dissipation field for these cases. In comparison to figure 
|l| we clearly see intermittent burst of turbulent activity in an otherwise smooth flow. 

The analysis is to some extent different from the homogeneous case. Firstly, 
most quantities like R\ and 77 loose their physical meaning. Secondly, the A^ plots 
become distorted by the laminar phases, so that the analysis is open to some greater 
degree of arbitrariness than in the homogeneous case. Table 2 gives the approximate 
values. 

From the spectra in Figure ^ it could be deduced that the small scale properties 
of the energy dissipation held are to some degree indepedent of the prerequisites of 
global uniformity and isotropy. We cannot state yet to which extent these condi- 
tions do not hold within a turbulent burst, but we accept the possibility that the 
multifractal aspects of the energy dissipaion held are in fact universal. 

Notice that the 'laminar' branch of the spectrum is not accessible due to the 
presence of a laminar background. Thus the right hand part of the spectrum col- 
lapses to a point at the maximum. The peak at unity indicates that the energy 
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Figure 4: The superposition of the f{a) curves of all locally isotropic experiments (Jl, J2, CFl, 
CF2) are in reasonable agreement. The solid line is the standard Jl, the errorbars given for CFl 
are representative for all experiments with 'poor' statistics. 
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5.1: The energy time series of experiment 5.2: The energy time series of experiment 
J3 in the off-centre region of a free jet. CVl in the wake close behind the cylinder. 
The width was chosen to match figure \i 
in terms of rj. Laminar phases are mixed 
with turbulent bursts which were trans- 
ported out of the central region of the jet. 
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flow 


free jet 


wake 




off-centre 


near a cylinder 


reference 


J3 


CVl 


position of probe 


60 diam. 
11.25 diam. off-centre 


4 diam. 


number of points 


1.25 X 10' 


1.25 X lO-" 


mean velocity 


0.43 m/s 


21.2 m/s 


sampling rate 


6 kHz 


200 kHz 


Kolmogorov 


0.17 mm 


0.035 mm 


rnicroscak^ r] 






Taylor niicroscalc A 


11 ;/ 


o3 // 


Taylor-Reynolds 


34 


720 


number Rx 






transitional range 


4 .. 25 r? 


18 .. 129 7] 


detector resolution 


6 r] 


28 T] 


permissible moments 


.. 4 


.. 4 



Table 2: Summary of experimental conditions and evaluation parameters for the in- 
homogeneous flows. 
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Figure 6: In the inhomogeneous case, the right branch of f(a) is not accessible. The solid 
line again is the standard Jl, errorbars are given for JS. Statistcs are quite poor for the jet 
measurement due to the great dilatation of the turbulent bursts. 
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dissipation is still space filling locally, despite its sparse appearance. 

6 Conclusion 

We have reported on a correspondence in small scale energy dissipation spectra in 
both locally isotropic and inhomogeneous turbulence. The results are plausible if 
one takes into consideration that firstly the analysis is restricted to small scales of 
less than IOO77, and secondly, by the observation that the inhomogeneous turbulence 
in the boundary layer of the free jet can be thought of as a composition of zones of 
turbulence transported out of the centre and inserts from the surrounding laminar 
flow □. Less accessible to intuition are our findings for the non-developed turbulent 
wake. 

To set our findings in context with previous work on small scale turbulence, 
we want to stress once more that our work focuses on the transitional regime, a 
range of length scales between the inertial range and the dissipative limit. Here 
we expect that the dissipation plays a role of increasing importance. Traditionally, 
scaling behaviour is investigated on larger scales in the inertial range, where the 
dissipation should be of minor influence. Some rpcent, work has shed new light 
on the turbulent cascade in the inertial range. ElEJOa. Our results suggest that 
in the neighbouring regime the statistics of the energy dissipation field become 
universally multifractal. The challenge remains to match the various regimes in 
which turbulence is described to finally arrive at a truly comprehensive view. 
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